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ABSTRACT 

We study the evolution of star formation activity of galaxies at 0.5 < 2; < 3.5 as a function of 
stellar mass, using very deep NIR data taken with Multi-Object Infrared Camera and Spectrograph 
(MOIRCS) on the Subaru telescope in the GOODS-North region. The NIR imaging data reach K ~ 
23-24 Vega magnitude and they allow us to construct a nearly stellar mass-limited sample down to 
~ 10 9 - 5-10 M even at z ~ 3. We estimated star formation rates (SFRs) of the sample with two 
indicators, namely, the Spitzer/MIPS 24/im flux and the rest-frame 2800 A luminosity. The SFR 
distribution at a fixed M star shifts to higher values with increasing redshift at 0.5 < z < 3.5. More 
massive galaxies show stronger evolution of SFR at z > 1. We found galaxies at 2.5 < z < 3.5 show a 
bimodality in their SSFR distribution, which can be divided into two populations by a constant SSFR 
of - 2 Gyr- 1 . Galaxies in the low-SSFR group have SSFRs of - 0.5-1.0 Gyr" 1 , while the high-SSFR 
population shows ~ 10 Gyr -1 . The cosmic SFRD is dominated by galaxies with M star = 10 10-11 M 
at0.5<z<3.5, while the contribution of massive galaxies with M sta r = 10 11-11,5 M shows a strong 
evolution at z > 1 and becomes significant at z ~ 3, especially in the case with the SFR based on 
MIPS 24 /im. In galaxies with M sta r = 10 10-11 - 5 M , those with a relatively narrow range of SSFR 
(< 1 dex) dominates the cosmic SFRD at 0.5 < z < 3.5. The SSFR of galaxies which dominate the 
SFRD systematically increases with redshift. At 2.5 < z < 3.5, the high-SSFR population, which 
is relatively small in number, dominates the SFRD. Major star formation in the universe at higher 
redshift seems to be associated with a more rapid growth of stellar mass of galaxies. 
Subject headings: galaxies: evolution — galaxies: high-redshift — infrared:galaxies 



1. INTRODUCTION 

Determining how galaxies build up their stellar mass 
is crucial for understanding galaxy formation and evo- 
lution. The star formation rate (SFR) and the stellar 
mass, which is approximately considered as the integral 
of the past SFR, are the most important physical proper- 
ties of galaxies to investigate their stellar mass assembly 
history. Many previous studies measured the comoving 
cosmic SFR density as a function of redshift, which pro- 
vides a global picture of star formation history in the 
universe, using the rest-frame UV luminosity (e.g., Lilly 
et al. 1996; Madau et al. 1996; Connolly et al. 1997; 
Madau et al. 1998; Cowie et al. 1999; Steidel et al. 1999; 
Schiminovich et al. 2005; Bouwens et al. 2007; Reddy et 
al. 2008), the mid-infrared (MIR) flux (e.g., Le Floc'h 
et al. 2005; Perez-Gonzalez et al. 2005; Babbedge et al. 
2006; Caputi et al. 2007; Magnelli et al. 2009), the optical 
nebular emission lines (Yan et al. 1999; Ly et al. 2007; 
Geach et al. 2008; Shim et al. 2009; Dale et al. 2010). 
These studies revealed that the cosmic SFR density in- 
creases by an order of magnitude from z ~ to z ~ 1, 
and then reaches a peak around z ~ 2 and decreases at 
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higher redshift (Hopkins 2004; Hopkins & Beacom 2006). 
On the other hand, recent near-infrared (NIR) surveys 
enable to measure the average stellar mass density of the 
universe and several studies based on deep NIR surveys 
found that a rapid evolution of the cosmic stellar mass 
density occurred at 1 < z < 3 (e.g., Dickinson et al. 
2003; Fontana et al. 2006; Kajisawa et al. 2009 and ref- 
erence therein). These results suggest that a significant 
fraction of the stellar mass in the present universe had 
been formed at z ~ 1-3. 

The next step to understand how star formation and 
stellar mass assembly proceeded in galaxies at 1 < z < 3 
is to investigate the SFR of these galaxies as a func- 
tion of stellar mass and redshift. Multi- wavelength data 
from surveys such as COMBO-17, GOODS, COSMOS, 
UKIDSS, and AEGIS, have been used to measure SFR 
and stellar mass of galaxies at 1 < z < 3 in many 
studies, where SFRs are estimated from various indica- 
tors, namely, the rest-frame UV luminosity (Juneau et 
al. 2005; Daddi et al. 2007; Feulner et al. 2007; Cowie 
& Barger 2008; Drory & Alvarez 2008; Mobasher et al. 
2009; Magdis et al. 2010), the Spitzer/MIPS 24 /im flux 
(Papovich et al. 2006; Elbaz et al. 2007; Noeske et al. 
2007; Zheng et al. 2007; Santini et al. 2009; Damen et 
al. 2009a; Damen et al. 2009b), the Ha emission line 
(Erb et al. 2006; Forster Schreiber et al. 2009), and the 
radio continuum emission (Pannella et al. 2009; Dunne 
et al. 2009). These studies found that the average SFR 
or SSFR of star- forming galaxies at a fixed stellar mass 
increases with redshift strongly. The rate of incline is 
similar in all stellar mass ranges investigated (Zheng et 
al. 2007; Damen et al. 2009a; Damen et al. 2009b). The 
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SSFR of star- forming galaxies at z ~ 2 is higher than 
that at z ~ 1 by factor of several and that at z ~ by 
factor of several tens (Daddi et al. 2007; Santini et al. 
2009; Pannella et al. 2009). In the studies at z ~ 2, how- 
ever, the sample selection is based on relatively shallow 
NIR/MIR data or deep but optical data, and therefore 
the sample becomes incomplete at relatively high stellar 
mass. For example, a depth of 23-24 AB magnitude in 
the NIR-MIR wavelength corresponds to a limiting stel- 
lar mass of - io 10 - 10 - 5 M at z ~ 2 (e.g., Dunne et al. 
2009). On the other hand, the selection in the optical 
band could preferentially pick up actively star-forming 
galaxies and could miss objects with relatively low SFR 
at low mass (e.g., Feulner et al. 2007; Kajisawa & Ya- 
mada 2006). If this is the case, the average SFR at 
low mass would be overestimated. In order to investi- 
gate the mass dependence of star formation activity of 
galaxies at z ~ 2 precisely, the selection based on more 
deep NIR/MIR data is desirable. 

On the other hand, the SFR as a function of stellar 
mass and redshift enables to estimate the contribution 
of galaxies in different stellar masses to the cosmic SFR 
density and its evolution. At z < 2, galaxies with stellar 
mass of ~ 10 10-11 M dominate the cosmic SFR density, 
and the decrease of the contribution of these galaxies 
from z ~ 2 to the present seems to result in the evolu- 
tion of the total SFR density of the universe (Cowie & 
Barger 2008; Mobasher et al. 2009; Santini et al. 2009). 
The evolution of the contribution to the cosmic SFR den- 
sity is also similar in all mass ranges at z < 1. Several 
studies suggest that the contribution of massive galax- 
ies with M st ar > 10 11 Mq increases with redshift more 
strongly than lower-mass galaxies at z > 1, although 
these massive galaxies do not dominate the total SFR 
density even at z ~ 2 (Juneau et al. 2005; Santini et al. 
2009). As mentioned above, however, only galaxies with 
relatively high stellar mass tend to be sampled in these 
studies, and therefore the evolution of the contribution of 
relatively low- mass galaxies is still unclear. Furthermore, 
the similar analysis for galaxies at z ~ 3 has not been 
performed due to the lack of sufficiently deep NIR/MIR 
data. Since the cosmic SFRD at z ~ 3 is estimated to be 
similarly high with that at z ~ 1-2 (Hopkins & Beacom 
2006), it is important to study how star formation occurs 
in galaxies with different stellar masses at z ~ 3. 

In this paper, we study the SFR as a function of stel- 
lar mass for galaxies at 0.5 < z < 3.5 in the GOODS- 
North field, using very deep NIR data from MOIRCS 
Deep Survey (MODS, Kajisawa et al. 2006; Ichikawa et 
al. 2007). The MODS data reach - 23-24 Vega magni- 
tude (~ 25-26 AB magnitude) in the K band, and they 
allow us to construct a stellar mass limited sample down 
to ~ io 9,5-10 M even at z ~ 3 and to investigate the 
distribution of the SFR at high redshift without biases 
for objects with very high SFR (approximately low stel- 
lar mass-to-light ratio) at low mass. We use the MODS 
data and publicly available multi- wavelength data of the 
GOODS survey to estimate redshift, SFR, and stellar 
mass of the sample. Section 2 describes the data set and 
the procedures of source detection and photometry. In 
Section 3, we construct the stellar mass limited sample. 
The methods for estimating the SFRs of the sample are 
presented in Section 4. We show the distribution of the 
SFR as a function of stellar mass and its evolution in 



Section 5, and discuss their implication in Section 6. A 
summary is provided in Section 7. 

We use a cosmology with Ho — 70 km s ^ Mpc ^, £l m — 
0.3 and = 0.7. The Vega-referred magnitude system 
is used throughout this paper, unless stated otherwise. 

2. OBSERVATIONAL DATA AND PHOTOMETRY 

We use the if-selected sample of the MODS in the 
GOODS-North region (Kajisawa et al. 2009, hereafter 
K09), which is based on our deep Ji^if s -bands imaging 
data taken with MOIRCS (Suzuki et al. 2008) on the 
Subaru telescope. Four MOIRCS pointings cover ~ 70% 
of the GOODS-North region (~ 103.3 arcmin 2 , hereafter 
referred as "wide" field) and the data reach J = 24.2, 
H = 23.1, K = 23.1 (5a, Vega magnitude). One of 
the four pointings is the ultra-deep field of the MODS 
(~ 28.2 arcmin 2 , hereafter "deep" field), where the data 
reach J = 25.1, H = 23.7, K = 24.1. A full description 
of the observations, reduction, and quality of the data is 
presented in a separate paper (Kajisawa et al. 2010). 

The source detection was performed in the if -band im- 
age using the SExtractor image analysis package (Bertin 
& Arnouts 1996). At first, we limited the samples to 
K < 23 and K < 24 for the wide and deep fields, 
where the detection completeness for point sources is 
more than 90%. Then we measured the optical-to- 
MIR SEDs of the sample objects, using the publicly 
available multi-wavelength data in the GOODS field, 
namely KPNO/MOSAIC (U band, Capak et al. 2004), 
HubbleSpaceTelescope/ Advanced Camera for Surveys 
(HST/ ACS; B, V, i, z bands, version 2.0 data; M. Gi- 
avalisco et al. 2010, in preparation; Giavalisco et al. 
2004) and Spitzer/IRAC (3.6/im, 4.5/im, 5.8/im, DR1 
and DR2; M. Dickinson et al. 2010, in preparation), as 
well as the MOIRCS J- and i^-bands images. Details 
of the multi-band aperture photometry are presented in 
K09. Following K09, we used objects which are detected 
above 2a level in more than two other bands in addition 
to the 5cr detection in the if -band, because it is difficult 
to estimate the photometric redshift and stellar mass of 
those detected only in one or two bands. The number 
of those excluded by this criterion is negligible (21/6402 
and 42/3203 for the wide and deep fields, respectively). 

3. SAMPLE SELECTION 

K09 estimated the redshift and stellar mass of the K- 
selected galaxies mentioned above and constructed a stel- 
lar mass-limited sample to study the evolution of the stel- 
lar mass function. We use the same stellar mass-limited 
sample in this study and describe briefly how the sample 
was constructed in the following. 

In order to estimate the photometric redshift and stel- 
lar mass, K09 performed spectral energy distribution 
(SED) fitting of the multi-band photometry described 
above (UBVizJHK, 3.6/im, 4.5/im, and 5.8/im) with 
population synthesis models. We here adopt the results 
with GALAXEV model (Bruzual & Chariot 2003). We 
used the model templates with exponentially decaying 
star formation histories with the decaying timescale r 
ranging between 0.1 and 20 Gyr and Calzetti extinction 
law (Calzetti et al. 2000) in the range of E(B-V) = 0.0- 
1.0. Metallicity is changed from 1/50 to 1 solar metal- 
licity. We assume Salpeter IMF (Salpeter 1955) with 
lower and upper mass limits of 0.1 and 100 M for easy 
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Table 1 

Number of galaxies in each redshift bin 



redshift 


Wide a 


Deep a 


combined a 


combined 




(103arcmin 2 K < 23) 


(28.2arcmin 2 K < 24) 




(MIPS 24/im-detected) 


0.5-1.0 


1618 (843) 


702 (306) 


1701 (843) 


626 


1.0-1.5 


1156 (348) 


481 (105) 


1241 (348) 


408 


1.5-2.5 


1008 (189) 


469 (75) 


1109 (189) 


417 


2.5-3.5 


307 (65) 


234 (47) 


370 (66) 


84 



a Number in the parenthesis indicates objects with spectroscopic 
redshift. 

comparison with the results in other studies. If we as- 
sume the Chabrier-like IMF (Chabrier 2003), the stellar 
mass is reduced by a factor of ~ 1.8. The model age is 
changed from 50 Myr to the age of the universe at the 
observed redshifts. The resulting photometric redshifts 
show a good agreement with spectroscopic redshifts (Fig- 
ure 1 in K09). If available, we adopted spectroscopic red- 
shifts from the literature (Cohen et al. 2000; Cohen 2001; 
Dawson et al. 2001; Wirth et al. 2004; Cowie et al. 2004; 
Treu et al. 2005; Chapman et al. 2005; Reddy et al. 2006; 
Barger et al. 2008; Yoshikawa et al. 2010), and performed 
the SED fitting fixing the redshift to each spectroscopic 
value for these galaxies. The estimated redshifts and 
stellar mass-to-luminosity (M/L) ratios from the best-fit 
templates were used to calculate the stellar mass. The 
uncertainty of the stellar mass is estimated taking into 
account of the photometric redshift error for those with- 
out spectroscopic redshift, and is discussed in detail in 
K09. 

K09 also investigated how the if -band magnitude limit 
affects the stellar mass distribution of the sample by us- 
ing the distribution of the rest-frame U — V color, which 
reflects the stellar M/L ratio well, as a function of stel- 
lar mass, and then determined the limiting stellar mass 
above which more than 90% of objects are expected to 
be brighter than the if -band magnitude limit as a func- 
tion of redshift (Figure 3 and 4 in K09). The limiting 
stellar mass is ~ 1 x 10 9 M (~ 3 x 10 8 M©) for the wide 
(deep) sample at z = 0.75, - 2 x 10 9 M (~ 8 x 10 8 
M ) at z = 1.25, - 5 x 10 9 M (~ 2 x 10 9 M ) at 
z = 2.0, - 1 x 10 10 M (~ 3 x 10 9 M©) at z = 3.0, re- 
spectively. We use only objects above this limiting mass 
in this study. Further details of determining the limiting 
stellar mass are given in K09. 

Table 1 lists the number of objects in each redshift bin. 
We use the same redshift bins as those used in K09. The 
bin width is sufficiently larger than the typical photo- 
metric errors, and each redshift bin includes a reasonable 
number of galaxies for studying the SFR as a function of 
stellar mass statistically. In the following analysis, we ba- 
sically use "combined" (wide + deep) sample except for 
the estimates of the comoving number density and SFR 
density (section 5.2 and 5.3), where we need to calculate 
V max separately for the wide and deep samples. 

4. ESTIMATE OF THE SFR 

In order to estimate the SFRs associated with the sam- 
ple galaxies, we use two SFR indicators, namely, rest- 
frame 2800 A luminosity and observed 24/im flux. The 
rest-frame UV luminosity probes the light from young 
massive stars, while the IR luminosity measures those 



photons absorbed and re-emitted by dust. We basically 
use a combination of the IR and (dust-uncorrected) UV 
luminosities for objects with Spitzer/MIPS 24/im detec- 
tion and adopt dust-corrected UV luminosity for those 
without the MIPS detection. For comparison, we also 
check the case where the dust-corrected UV luminosity 
is used for all galaxies. 

4.1. SFR based on MIPS 2^m flux 

Public MIPS 24/im data for the GOODS-North 
(DR1+, M. Dickinson et al., in preparation) was used for 
the MIR photometry. Since this image is very deep and 
its PSF size is relatively large (FWHM ~ 5.4 arcsec), the 
source confusion occurs for many objects. Following Le 
Floc'h et al. (2005), we use the IRAF/DAOPHOT pack- 
age (Stetson 1987) in order to deal with these blended 
sources properly. The DAOPHOT software is based on 
the PSF fitting technique and allows us to fit blended 
sources in a crowded region simultaneously and then de- 
rive the flux of each object from the scaled fitted PSF. 
An empirical PSF was constructed from bright isolated 
point sources. We used the positions of the sample galax- 
ies on the high-resolution MOIRCS if s -band image as a 
prior for the centers of the fitted PSFs in the photome- 
try. The background noise was estimated by measuring 
sky fluxes at random positions on the image. For sources 
whose residual of the PSF fitting was significantly larger 
than the background noise, we added the fitting residual 
to the photometric error. We set a detection threshold of 
5<j, which corresponds to ~ 20/iJy in most cases where 
the residual of the fitting is negligible. 

In order to convert the MIPS 24/im flux to the total IR 
luminosity (8-1000 /im), we used SED templates of star- 
forming galaxies provided by Dale & Helou (2002). The 
model templates cover a wide range of SED shapes, al- 
lowing for different heating levels of the interstellar dust 
by the radiation field. Following Wuyts et al. (2008) and 
Damen et al. (2009a), we calculated the total IR luminos- 
ity for each object using the all templates with the range 
from a = 1 to a = 2.5. a is defined by dM(U) ~ U~ a dU, 
where M(U) represents the dust mass headed by an in- 
tensity U of the radiation field (Dale & Helou 2002). 
Then we averaged the calculated values for the different 
a in a logarithmic scale and adopted as a best estimate 
for the IR luminosity. The variation from the value for 
a = 1 to that for a = 2.5 was taken into account in the 
computation of the error of the total IR luminosity. For 
objects without spectroscopic redshift, we also estimated 
the effect of the uncertainty of the photometric redshift 
by calculating the IR luminosity over the possible ranges 
of redshift, and added it to the error. 

We converted the derived total IR luminosity to a 
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Figure 1. Comparison between SFRir + uv an d SFR correcte d_uv f° r galaxies detected in the MIPS 24/im image for each redshift bin. 
Squares represent X-ray sources detected in the Chandra hard-band image. Typical errors of the SFRjr+uv and SFR correcte d— tjv are 
shown at the bottom right corner of each panel. 



SFRir using the Kennicutt (1998) 's relation 

SFR lR (M yr' 1 ) = 4.5 x 1(T 44 L lR (erg s' 1 ). (1) 

Salpeter IMF with lower and upper mass limits of 0.1 
and 100 M is assumed. The limiting flux of ~ 20/iJy 
mentioned above corresponds to ~ 7 M /yr at z = 1, 
- 40 M /yr at z = 2, and - 300 M /yr at z = 3. The 
number of objects detected in the MIPS 24 /im image in 
each redshift bin is 626 at 0.5 < z < 1.0, 408 at 1.0 < 
z < 1.5, 417 at 1.5 < z < 2.5, and 84 at 2.5 < z < 3.5, 
respectively. 

4.2. SFR based on rest-frame 2800 A luminosity 

The rest-frame 2800 A luminosity for each object was 
calculated from the same best-fit SED template as used 
in the estimate of the photometric redshift and stellar 
mass. Since our multi-band photometric data cover com- 
pletely the rest-frame 2800 A for 0.5 < z < 3.5, the esti- 
mated rest-frame 2800 A luminosity is a robust quantity, 
especially in the case where a galaxy is detected in both 
bands which straddle the rest-frame 2800 A. We calcu- 
lated the confidence interval for it in the SED fitting 
procedure. In order to derive the dust-corrected UV lu- 
minosity, we also used E(B-V) of the same best-fit SED 
template. Since E(B-V) is often degenerated with other 



parameters of stellar population such as age and metal- 
licity in the SED fitting, the dust-corrected UV luminos- 
ity tends to have larger errors than the dust-uncorrected 
one. Note that our lower limit of the model age of 50 Myr 
in the SED fitting affects the estimated E(B-V), thus the 
dust-corrected UV luminosity and SFR. We will discuss 
the effect of the age limit on SFR in Section 5.1. 

We converted the dust-uncorrected and dust-corrected 
rest-frame 2800 A luminosities into SFRs using the cali- 
bration by Kennicutt (1998) 

SFRw (M yr' 1 ) = 1.4 x 10~ 28 L v (erg s~ x Hz' 1 ). 

(2) 

The SFR estimated from the dust-uncorrected UV lumi- 
nosity accounts for the unobscured light from young stars 
and is added to the SFRir estimated from the 24/im flux 
for objects detected in the MIPS 24/im image. We re- 
fer to the combined value as SFRir + tjv in the following. 
The SFR estimated from the dust-corrected UV lumi- 
nosity (hereafter SFR correcte d-uv) can be calculated for 
all the sample. 

4.3. Comparison between 5Fi?iR + uv and 

SFi^corrected-UV 

Figure 1 shows a comparison between SFRir + tjv and 
SFR correc ted— uv for galaxies detected in the MIPS 24/im 
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Figure 2. SFR vs. M s tar for our stellar mass limited sample for each redshift bin. Squares show galaxies with SFRir + uv derived from the 
24 /mi fluxes, while circles represent objects undetected in the MIPS 24 fim image, with SFR correctec i_uv derived from the dust-corrected 
rest-frame UV luminosity. Median errors of the SFR and stellar mass at each mass are shown as the crosses at the bottom of each panel. 
Dotted lines represent SFR/M sta r = 1/50 Myr -1 (see text for details). Double-dotted dashed lines show SFR/M sta r = 0.1, 1, 10 Gyr -1 . 
Blue dashed lines show the median SFR at each stellar mass. Green dotted-dashed lines represent the results of the linear fitting for 
star-forming galaxies. In the panel of the 2.5 < z < 3.5 bin, the lower line shows the result for galaxies with log(SFR/M s tar) > -0.5, and 
upper one shows that for galaxies with log(SFR/M s tar) > 0.25. 



image. Given the uncertainty, these two estimates agree 
well with each other especially in low-redshift bins. At 
low SFR, however, SFRir + tjv tends to be higher than 
SFR correcte d-uv, and the systematic offset seems to be- 
come larger with redshift. The mean and standard de- 
viation of log(SFR m+V v/SFR corrected ^uv) for all ob- 
jects with the MIPS 24/im detection are 0.19 ± 0.39 at 
0.5 < z < 1.0, 0.25 ± 0.37 at 1.0 < z < 1.5, 0.37 ± 0.49 
at 1.5 < z < 2.5, and 0.63 ± 0.64 at 2.5 < z < 3.5. Franx 
et al. (2008) reported a similar trend in the comparison 
between the SFRs derived from the MIPS 24/im flux and 
the SED fitting. 

In the conversion from the 24/im flux into the total IR 
luminosity, we don't take into account the luminosity- 
dependent shapes of the IR SED found in local galaxies 
(e.g., Chary & Elbaz 2001), following Wuyts et al. (2008) 
and Damen et al. (2009a). Therefore SFRir + tjv could be 
overestimated for galaxies with low IR luminosity (SFR) 
and underestimated for those with high IR luminosity, if 
the luminosity-dust temperature (IR color) relation for 
nearby galaxies holds for galaxies at higher redshift. On 
the other hand, there are several studies that suggest 



that high-z active star-forming galaxies do not form such 
a relation and galaxies with Lir > 1O 12 L (ULIRGs) at 
high redshift have cooler average dust temperatures than 
those in the local universe (e.g., Symeonidis et al. 2009; 
Muzzin et al. 2010). Furthermore, it is reported that 
bright 24 /im sources at z ~ 2 have higher polycyclic 
aromatic hydrocarbon (PAH) luminosity relative to the 
total IR luminosity than the local ULIRGs with similar 
IR luminosity (Huang et al. 2009; Murphy et al. 2009). 
Since the MIPS 24 /im band samples the PAH emissions 
in 6-9 /im region at z > 1.5, higher Lpah/Lir ratios than 
those for the local ULIRGs could lead to a overestima- 
tion of SFR for these high-z galaxies if one converts 24 
/im fluxes into the SFRs using the local luminosity-dust 
temperature relation. 

It should be noted that the MIR flux may also 
be powered partly by AGN (e.g., Nardini et al. 
2008). We showed objects detected in the Chan- 
dra hard-band image of the CDF-North (Alexander et 
al. 2003) as squares in Figure 1. Those objects de- 
tected in the hard-band of Chandra have slightly higher 
SFRiR + uv/SFR correct ed-uv ratios; the mean and stan- 
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Figure 3. The same as Figure 2 except that SFR, 
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dard deviation of log(5Fi?iR + uv/S , ^ 1 -Rcorrected-uv) for 
those objects are 0.31 ±0.38 at 0.5 < z < 1.0, 0.37 ±0.35 
at 1.0 < z < 1.5, 0.54 ± 0.32 at 1.5 < z < 2.5, and 
0.80 ± 0.66 at 2.5 < z < 3.5. SFRt R+uv for those 
sources could be overestimated due to the AGN con- 
tribution, although SFRir + uv tends to be higher than 
SFR correcte d-uv especially at high redshift even if we ex- 
clude the Chandra hard-band sources. 

On the other hand, the dust-corrected UV luminosity 
might be underestimated for heavily obscured galaxies 
such as ULIRGs at low and high redshifts (e.g., Goldader 
et al. 2002; Papovich et al. 2006; Reddy et al. 2010). If 
a galaxy has star-forming regions from which one can 
detect no UV light at all due to the heavy dust obscura- 
tion, only UV light from relatively less-obscured regions 
contributes the observed SED, which could result in the 
underestimate of the dust extinction. 

Since the calibration by Kennicutt (1998) we used for 
the conversion from the rest-frame 2800 A luminosity 
into SFR is derived assuming that SFR has remained 
constant over a timescale of ~ 100 Myr and that the 
luminosity is dominated by young O and B stars, SFR 
could be overestimated for relatively quiescent galaxies. 
For such galaxies, older stellar population than O and 
B stars could contribute significantly to the rest-frame 
2800 A luminosity (e.g., Dahlen et al. 2007). Therefore 



SFR correcte d-uv could be systematically overestimated 
for galaxies with low SFRs relative to their stellar mass. 

Keeping in mind the possible systematic errors, we 
show both results with SFRir + tjv for objects detected 
in the MIPS 24/im image and SFR correcte d-uv for those 
without the MIPS detection, and with SFR correct ed-uv 
for all the sample in the following section. 

5. RESULTS 

5.1. SFR as a function of Stellar Mass at 0.5 < z < 3.5 

Figure 2 and 3 show the SFR as a function of M sta r 
for the stellar mass selected sample in the four redshift 
bins. Figure 2 represents the case where SFRir + tjv 
is used if a galaxy is detected in the MIPS 24/im im- 
age and SFR correc ted-uv is used for those without the 
24/im detection, and Figure 3 represents the case where 
SFR correc ted-uv is used for all the sample. A dot- 
ted line in each redshift panel represents SFR/M sta r = 
1/50 Myr -1 . In the calculation of the SFR correc ted-uv 7 
both the rest-frame 2800 A luminosity and the E(B-V) 
were estimated from the SED fitting. The stellar mass 
is also based on the same best-fit SED model. Since 
SFR/M 

star of the model SED cannot exceed the inverse 
of the age for the monotonically decaying SFR mod- 
els (SFR/M 

star roughly equals the inverse of the age 
for the constant SFR models if one ignores the stellar 
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Figure 4. SFR distribution for galaxies in each stellar mass and redshift bin. SFRir+tjv is used if a galaxy is detected in the MIPS 
24/mi image and SFR correcte d-UV is used for objects without the 24/i,m detection. Dashed and dotted-dashed lines show the median and 
average SFRs for each mass and redshift bin. 



mass loss due to the super nova or stellar wind), the 
lower limit of the model age leads to the upper limit for 
the SFR/M 

star- Therefore, the lower limit of age > 50 
Myr in the SED fitting mentioned above causes the up- 
per limit of SFR corrected _ uv /M star < 1/50 Myr" 1 (= 20 
Gyr -1 ). On the other hand, the SFRir + tjv is not af- 
fected by this limit, because the SFR based on the 24 
/am flux are independent of the result of the SED fitting 
for the stellar emission. Since there are only a few ob- 
jects whose SFRir + uv significantly exceeds this limit in 
high-redshift bins in Figure 2, we believe that the lower 
limit of the model age does not significantly bias the dis- 
tribution of the SFR correcte d-uv as a function of M sta r- 

In Figure 2 and 3, we plot the median and average of 
SFR at each mass as a dashed and dotted-dashed lines, 
respectively, to examine in detail the stellar mass depen- 
dence of the SFR. At M star < 10 10 5 M , more massive 



galaxies tend to have higher SFRs in all redshift bins, 
although there is a substantial scatter at each mass. The 
median and average SFRs of these galaxies increase with 
M s tar- On the other hand, the trend becomes less clear 
at M s tar > 10 10 - 5 M ;the median SFR does not signif- 
icantly increase with M s tan especially in the case with 
SFR correc ted-uv for all the sample (Figure 3). There are 
massive galaxies with relatively low SFRs similar with 
low-mass galaxies, while galaxies with higher SFRs which 
follow the trend at lower M s tar also exist. 

We note that the SFR distribution at a fixed M s tar 
shifts to higher values with increasing redshift at least at 
z < 2.5 in Figure 2 and 3. Figure 4 and 5 show the SFR 
distribution for each stellar mass and redshift bin. At 
z < 2.5, the shift of the SFR distribution with redshift is 
seen in all mass bins in the both figures, although there 
are galaxies with relatively wide range of SFR in each 
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Figure 5. The same as Figure 4 except that SFR correcte d_uv is used for all the sample. 



mass and redshift bin. Furthermore, it is seen that more 
massive galaxies show stronger evolution of the SFR in 
the both cases with different SFR indicators. 

At 2.5 < z < 3.5, galaxies with M star = 10 9 5 -10 n M 
show a bimodality in the SFR, and the SFR distribution 
does not concentrate around the median or average value 
(Figure 4 and 5). The SFRs of the both high-SFR and 
low-SFR groups of the bimodality seem to increase with 
M s tar- Since the bimodality is also seen in the case where 
the SFR correcte d-uv is used for all the sample, it is not 
due to the effect of the detection limit of the MIPS 24 jam 
data. While the median SFRs of galaxies at 2.5 < z < 3.5 
are similar with those at 1.5 < z < 2.5 in all stellar mass 
bins, the average SFRs of these galaxies tend to be higher 
than those atl.5<z<2.5 due to galaxies with very high 
SFRs, especially at high mass. 

5.2. SSFR vs M star 



In Figure 6 and 7, we show SSFR vs M s tar for our 
sample in the both cases with the different SFR indica- 
tors. At Mgtar ^ 10 10 5 Mq, the median and average 
of the SSFR are nearly independent of M s tar at least at 
z > 1, although the median and average values slightly 
decrease with M s tar in the lowest redshift bin. This indi- 
cates that the relation between SFR and M s tar for galax- 
ies with Mgtar ^ 10 10 ' 5 M has a logarithmic slope of ~ 
1 (slightly less than 1 for the lowest redshift bin). On 
the other hand, the median and average SSFRs decrease 
with Mgtar at Mgtar ^ 10 10 ' 5 M , except for the average 
SSFR in the 2.5 < z < 3.5 bin. 

In order to investigate the evolution of the SSFR, we 
show a distribution of SSFR for the cases with the differ- 
ent SFR indicators in Figure 8. We use the V max method 
to calculate the comoving number density. K09 calcu- 
lated the V m ax for the if -band magnitude limit (K = 23 
for the wide field or K = 24 for the deep field) by using 
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Figure 6. Specific star formation rate (SFR/M s tar) vs M s tar for our stellar mass limited sample for each redshift bin. SFRir+tjv is used 
if a galaxy is detected in the MIPS 24/mi image and SFR correcte d-UV is used for objects without the 24/xm detection. Symbols are the 
same as Figure 2. 



the best-fit model template for each object (see Section 
3.3 in K09 for details). In the figure, we use only galax- 
ies with M star > 10 9 ' 5 M for a fair comparison among 
the different redshift bins. We also show galaxies with 
M star = 10 9 - 5 - 10 - 5 M and with M star > 10 10 5 M sepa- 
rately. It is seen that the distribution of the SSFR shifts 
to higher values with redshift, while the comoving num- 
ber density of galaxies decreases. 

At2.5<z<3.5, the bimodality in the SFR mentioned 
above is also seen in the SSFR distribution, as a more 
simple form. There are a group of galaxies with high 
SSFRs of SFR/M s t ar - 10 Gyr" 1 and that of galaxies 
with relatively low SSFRs of SFR/M sta r ~ 0.6 Gyr" 1 . 
Thus, a constant SSFR of SFR/M sta r ~ 2 Gyr" 1 can 
divide between the two populations independent of M s tar 
in the both cases with the different SFR indicators. 

5.3. Contribution to the Cosmic SFR Density 

Here we examine the contribution to the cosmic SFR 
density (hereafter, SFRD) of galaxies with different stel- 
lar masses. We use the same V ma x as used in Section 5.2 
to estimate the SFRD. At first, we show the integrated 
SFRD for our sample and compare it with a compilation 
by Hopkins (2004) in Figure 9. The integrated SFRD 
was calculated from the sample limited by the stellar 



mass which depends on redshift (Section 3) and was not 
extrapolated down to lower mass. Therefore the SFRD 
at higher redshift in the figure was integrated only down 
to the higher stellar mass limit. Since the contribution 
of galaxies with M s t a r < 10 10 M seems to be negli- 
gible even at high redshift as we show below, however, 
the effect of the systematically higher mass limit at high 
redshift is expected to be small. 

In Figure 9, we plot the both results with the different 
SFR indicators and also show the case with SFRir + tjv 
where Chandra hard-band sources are excluded. The 
integrated SFRDs for our sample increase with redshift 
up to z ~ 2 and are constant or slightly decrease between 
z ~ 2 and z ~ 3. They range within the compilation by 
Hopkins (2004) over 0.5 < z < 3.5, except for a lower 
SFRD at z > 1.5 in the case with SFR correcte d-uv for all 
the sample. S F R correcte d-uv f° r galaxies at high redshift 
might be underestimated due to the heavy obscuration 
(Section 4.3). 

Figure 10 shows the contribution to the cosmic SFRD 
of galaxies in the different ranges of stellar mass. In 
the all redshift range, galaxies with M s tar = 10 10-11 M 
dominate the cosmic SFRD; their contribution is about 
1/2 - 2/3 of the total SFRD. The increase of the cosmic 
SFRD with redshift up to z ~ 2 seems to reflect mainly 
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Figure 7. The same as Figure 6 except that SFR, 
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is used for all the sample. 



the evolution of the contribution of galaxies in this mass 
range. 

In the case with SFRtr + tjv, the contribution of mas- 
sive galaxies with M sta r = 10 11-11 ' 5 M shows stronger 
evolution at z > 1.5 than that of galaxies with M sta r = 
-LQio-ii jyj-^ anc j continues to increase up to z ~ 3. 
These massive galaxies significantly contribute the cos- 
mic SFRD at 2.5 < z < 3.5, while the contribution of 
these galaxies is relatively small at z < 2. In the case 
where the Chandra hard-band sources are excluded, the 
contribution of these massive galaxies is slightly lower 
and its evolution is milder, although the contribution 
becomes significant at z ~ 3 even in this case. This 
is because the X-ray AGN is preferentially associated 
with massive galaxies (e.g., Best et al. 2005) and be- 
cause the fraction of objects detected at hard X-ray in 
massive galaxies becomes very high at z > 2 (Yamada et 
al. 2009; Brusa et al. 2009). While the SFRt R+uv of the 
Chandra hard-band sources could be overestimated due 
to the contribution of AGN, the exclusion of these galax- 
ies might cause the underestimate of the contribution of 
massive galaxies to the SFRD if a significant fraction of 
the IR luminosity of these galaxies arises from star for- 
mation activities. In the case where SFR correcte d-uv is 
used for all the sample, the contribution of these galax- 
ies with M s tar = 10 11-11 - 5 M has a large uncertainty 



even at low redshift and its evolution is less clear. This 
is partly because the large error of the SFR for heavily 
obscured galaxies with high SFR. In particular, the con- 
tribution at 1.0 < z < 1.5 is strongly affected one such 
object with SFR corrected _ uv > 1000 M /yr (Figure 3). 

The contribution of galaxies with M sta r < 10 10 M 
seems to be small at any redshift, although the complete- 
ness limit for our sample becomes ~ iq 9,5-10 y[ a t z ~ 
3. The contribution of galaxies with M star = 10 9 ^ -10 M 
is - 15% of the total SFRD at 0.5 < z < 3.5 and evolves 
similarly with that for galaxies with M s tar = 10 10-11 M . 
Galaxies with M s tar > 10 11,5 M have a negligible contri- 
bution, although there are only a few such very massive 
galaxies in each redshift bin and the uncertainty of the 
contribution is very large. 

6. DISCUSSION 

6.1. Comparison with Other Studies at z < 2.5 

We have investigated the distribution of SFR as a func- 
tion of stellar mass for galaxies at 0.5 < z < 3.5, using 
the stellar mass limited sample based on the very deep 
NIR imaging data of the MODS. We here compare our 
results in the MODS field with previous studies mainly 
at z < 2.5. 

Several studies examined the relation between SFR 
and stellar mass of star- forming galaxies, and found that 
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Figure 8. Distribution of the specific star formation rate 
(SFR/Mstar) for galaxies with M sta r > 10 95 M in each redshift 
bin. Left panels show the case where SFRir+tjv is used for galax- 
ies detected at 24 fim and SFR correcte d_uv for the other objects, 
while right panels show the case where SFR corrected _uv is used 
for all the sample. Long dashed line (Magenta) represents galaxies 
with M s tar > 10 10 - 5 Mq, and short dashed line (blue) shows those 
with M s tar = 10 9 - 5-10 - 5 Mq. The insets represent the zoom-up of 
the plots for galaxies with M s tar > 10 10 5 Mq. 



more massive galaxies tend to have higher SFRs at z ~ 1 
(Noeske et al. 2007; Elbaz et al. 2007) and at z ~ 2 (San- 
tini et al. 2009; Daddi et al. 2007; Pannella et al. 2009; 
Dunne et al. 2009). At M star < 10 10 5 M , we found 
a similar trend in our sample at 0.5 < z < 2.5, and 
confirmed the same trend at 2.5 < z < 3.5 although 
the SFR distribution shows the bimodality and does not 
concentrate around the average value. At M star > 10 10 5 
Mq, there are massive galaxies with similar SFRs with 
lower-mass galaxies in our stellar mass-selected sample, 
while galaxies with higher SFRs which follow the trend 
at lower M sta r also exist. If we excluded galaxies with 
relatively low SFRs as a quiescent population, we could 
consider the relation between SFR and stellar mass con- 
tinues to high mass for star- forming galaxies, which has 
been reported in other studies. 
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mass limited sample for the cases with the different SFR indicators. 
The case with SFRjr+uv wnere the Chandra hard-band sources 
are excluded is also shown (the result with SFR correc t e d— uv does 
not change significantly if we exclude the Chandra hard-band 
sources). Small open circles show the compilation from Hopkins 
(2004). 



Our result that the SFR (SSFR) distribution at a fixed 
M s tar shifts to higher values with increasing redshift is 
also consistent with these previous studies. For exam- 
ple, Pannella et al. (2009) estimated an average SFR as 
a function of M sta r for star-forming BzK galaxies with 
K < 21.2 at z ~ 2 in the COSMOS field by performing 
a stacking analysis of 1.4 GHz radio continuum. They 
pointed out that the SSFR of these star- forming galaxies 
show a evolution of a factor of ~ 4 between z = 1.4 and 
z = 2.3. Using the stellar mass-selected sample whose 
limiting mass reaches ~ 10 9 ' 5 -10 10 M even at z ~ 3, 
we found that more massive galaxies show stronger evo- 
lution of the SFR at z > 1. 

We found the median and average SSFRs of galaxies 
is nearly independent of M sta r at M sta r ^ 10 10 5 M at 
z > 1, while they slightly decrease with M sta r at 0.5 < 
z < 1.0. This indicates that the logarithmic slope of 
the relation between SFR and M star is n ~ 1 (n < 1 
for 0.5 < z < 1.0). Other studies reported the similar 
slope of the SFR-M star relation, although there is some 
variance among the studies in the observed slope and 
scatter of the relation (Daddi et al. 2007; Santini et al. 
2009; Pannella et al. 2009; Elbaz et al. 2007; Noeske et 
al. 2007). The variance might be explained by differences 
in the sample selection or the SFR indicators among the 
studies. 

The SSFR independent of M star at M star < 10 10 5 M 
means that the stellar mass growth by star formation 
activities in each galaxy is similar among galaxies with 
different masses. Using the same sample as this study, 
K09 investigated the evolution of the galaxy stellar mass 
function (SMF), and found that the low-mass slope of the 
SMF becomes steeper with redshift at 1 < z < 3. The 
nearly constant mass growth rate by star formation ex- 
pected from the slope of n ~ 1 implies that star formation 
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Figure 10. Contributions of galaxies in different stellar mass ranges to the cosmic SFR density as a function of redshift. Upper left and 
right panels show the cases with the different SFR indicators. Bottom left panel represents the case with SFRjr+uv where the Chandra 
hard-band sources are excluded. Open symbols mean that the incompleteness due to the stellar mass limit could be significant for the 
redshift bin. 



in each galaxy does not significantly change the low-mass 
slope with time. If such mass-independent SSFR contin- 
ues down to lower mass than the stellar mass limit of our 
sample in the high redshift bins, additional mechanisms 
such as the hierarchical merging, which destroys small 
galaxies and builds more massive ones, might be needed 
in order to reproduce the evolution of the low-mass slope 
(K09). 

6.2. Bimodality in SSFR of star-forming galaxies at 
2.5 < z < 3.5 

We found the bimodality in the SSFR distribution of 
galaxies at 2.5 < z < 3.5, which can be divided into 
two populations by a constant SSFR of ~ 2 Gyr -1 . One 
population consists of galaxies with SFR/M st ar ~ 0.5—1 
Gyr -1 , and those with SFR/M star - 10 Gyr -1 belong to 
the other population. For the SSFR of galaxies at z > 
2.5, the SFR and M sta r of Lyman Break Galaxies (LBGs) 



have been investigated by several studies. Magdis et al. 
(2010) reported that LBGs detected in the IRAC 3.6 
fim and 4.5 fim bands at z ~ 3 show the SFR-M sta r 
relation with a slope of n ~ 0.9 and a considerable scatter, 
and that their median SSFR of ~ 4.5 Gyr -1 is higher 
than those at lower redshift s, while several observational 
studies suggest that the SFR-M star relation for Lyman 
Break Galaxies at 4 < z < 7 is similar with that at 
z~2 (Daddi et al. 2009; Stark et al. 2009; Gonzalez et 
al. 2010). The range of the SSFR of LBGs in Magdis 
et al. (2010) is similar with galaxies at 2.5 < z < 3.5 in 
our sample, although most of LBGs in their sample have 
M star > 10 10 - 5 M probably due to the IRAC selection 
and the number of galaxies with low SSFR is relatively 
small. 

Dutton et al. (2010) constructed a semi-analytic model 
for disk galaxies and suggested that the evolution of the 
normalization of the SFR-M sta r relation for these disk 
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Figure 11. Rest-frame broad-band SEDs of star-forming galaxies at 2.5 < z < 3.5. Left and middle panels show those with 
log(SFRi R+uv /M st ar) > 0.25 (the high-SSFR population) and those with log(SFRi R+uv /M st ar) = -0.5 - 0.25 (the low-SSFR popu- 
lation). Errorbars represent the la photometric uncertainties and arrows indicate the 2a upper limits. Solid line in each panel shows the 
best-fit SED model for all the data points. Average stellar age, color excess, and metallicity of the best-fit model are also shown in the pan- 
els. Right panel shows a comparison between these average SEDs for the high- and low-SSFR populations. If we use the SFR correctec i_uv 
for all the sample, the best-fit SEDs for the high- and low-SSFR populations are almost the same as seen in this figure. 



galaxies at z < 2 closely follows the evolution in the gas 
(and dark matter) accretion rate. They predicted that 
the normalization continues to increase with redshift up 
to z > 2, if the SFR follows the gas accretion rate at such 
high redshift. The SFR of the high-SSFR population at 
2.5 < z < 3.5 may be related with the evolution of the 
gas accretion rate at such high redshift, although the 
SSFR of ~ 10 Gyr -1 at z ~ 3 is higher than the model 
prediction by Dutton et al. (2010). 

Recently, Daddi et al. (2010) reported a bimodality in 
the SFR relative to the gas mass for star-forming galax- 
ies at low and high redshifts and proposed that the star 
formation efficiency would be different between disk and 
starburst galaxies because different modes of star forma- 
tion work in the two populations. The high-SSFR pop- 
ulation at 2.5 < z < 3.5 in our sample might correspond 
to such starburst galaxies. For example, the starburst 
sequence consists of sub-millimeter galaxies (SMGs) at 
high redshift as well as local ULIRGs in Daddi et al. 
(2010) and several studies suggest that high-z SMGs tend 
to show relatively high SSFRs (e.g., Takagi et al. 2008; 
Daddi et al. 2009). There are 23 SCUBA sources in the 
MODS field (Pope et al. 2006), and all but one (GN10, 
Wang et al. 2007; Wang et al. 2009) have counterparts 
in the if -selected catalog of the MODS. Of 22 sources 
detected in the if -band image, 6 SCUBA sources lie at 
2.5 < z < 3.5 (3 with spectroscopic redshift and 3 only 
with photometric redshift), all but one hard X-ray source 
(GN22, Laird et al. 2010) belong to the high-SSFR pop- 
ulation with SFR/M s tar ~ 10 Gyr -1 . Figure 11 shows 
the rest-frame broad-band SEDs of star-forming galaxies 
at 2.5 < z < 3.5 on the high- and low-SSFR popula- 
tions. We defined the high-SSFR population as galaxies 
with lo^SFR/MstarlGyr- 1 ]) > 0.25 and the low-SSFR 
population as those with log(SFR/M star ) = -0.5-0.25, 
respectively. The high-SSFR population shows weaker 
B aimer break feature and redder color at long wave- 
length, which indicates that these galaxies are younger 
and dustier, while the overall slope of SEDs and scatters 
among the objects are similar between the two popula- 
tions. Those objects in the high-SSFR population might 



have higher star formation efficiency than ordinary star- 
forming (disk) galaxies. 

6.3. Evolution of Cosmic SFRD and Stellar Mass 
Assembly of Galaxies 

In Section 5.3, we saw that galaxies with M sta r = 
10 10 - 11 M dominate the cosmic SFRD at 0.5 < z < 3.5, 
and the contribution of galaxies with M sta r = 10 11-11 ' 5 
Mq increases with redshift even at z > 2 and becomes 
significant at z ~ 3. On the other hand, K09 inves- 
tigated the contribution of galaxies in different mass 
ranges to the cosmic stellar mass density, and found 
that galaxies with M star - iqIO.5-11.5 Mq ( 10 io-ii.5 
Mq at z ~ 3) dominate the cosmic stellar mass den- 
sity (Figure 17 in K09). Considering the nearly con- 
stant SSFR among (star- forming) galaxies with different 
stellar masses discussed in Section 6.1, the result that 
galaxies with io 10-11 - 5 Mq dominate the cosmic SFRD 
is qualitatively consistent with the dominant contribu- 
tion of these galaxies in the stellar mass density. Major 
star formation in the universe at 1 < z < 3 seems to 
occur in galaxies with stellar mass which dominates the 
cosmic stellar mass density or those with slightly lower 
stellar mass. 

The cosmic SFRD increases with redshift from z ~ 
up to z ~ 2 and is roughly constant or slightly decreases 
between z ~ 2 and z ~ 3 (e.g., Hopkins & Beacom 2006). 
The evolution of the cosmic SFRD roughly follows that 
of galaxies with M star = 10 10-11 Mq, which dominate 
the SFRD (Figure 9 and 10). In Figure 4 and 5, the 
average SFR per galaxy for those with M star = 10 10-11 
M increases with redshift up to z ~ 3, but the evolu- 
tion between z ~ 2 and z ~ 3 seems to be weaker than 
that at lower redshift. On the other hand, K09 inves- 
tigated the evolution of the number density of galaxies 
in different mass ranges, and showed that the number 
density of these galaxies decreases gradually with red- 
shift (Figure 16 in K09). At z < 2, the increase in the 
average SFR per galaxy with redshift overwhelms the 
decrease in the number density, and therefore the con- 
tribution of these galaxies to the total SFRD increases 
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Figure 12. Contributions of galaxies in different SSFR ranges 
to the cosmic SFR density in each redshift bin for the case 
where SFRir + uv is used for galaxies detected at 24 fim and 
SFR correcte d_ uv for the other objects. Only galaxies with M s tar > 
10 9-5 Mq are used for all the redshift bins. Thick solid line shows 
galaxies with M s t ar > 10 9-5 Mq and thin solid line represents 
galaxies with M sta r = 10 10-115 M . Dotted (magenta), long 
dashed (green), and short dashed (blue) lines represent galaxies 
with M s tar = 10 11 - 115 , 10 10 - 5 - 11 , and lO 10 " 10 - 5 M , respec- 
tively. 

with redshift. At 2 < z < 3, the evolution of the SFR 
per galaxy becomes weaker, while the number density 
continues to decrease even at such high redshift. Then 
these two roughly balances, and the SFRD is nearly con- 
stant or slightly decreases with redshift. For galaxies 
with M star = 10 11-11 - 5 M , the average SFR evolves 
more strongly between z ~ 2 and z ~ 3 especially in the 
case with SFRtr+tjv- The strong evolution of the SFR 
per galaxy for these massive galaxies continues to over- 
whelm the decrease in the number density up to z ~ 3, 
although the number density of these galaxies also con- 
tinues to decrease to z ~ 3. As a result, the contribution 
of these galaxies with M star = lO 11-11 - 5 M to the SFRD 
increases with redshift even at z > 2. 

In other words, the number density of galaxies is rela- 
tively small, but the average SFR per galaxy is high at 
z ~ 3. As time goes on, the SFR in each galaxy de- 
creases, while the number density of galaxies at a fixed 
M s tar increases due to the stellar mass growth in galax- 
ies with a wide range of masses. In particular, the SFR 
in massive galaxies with M s tar = 10 11-11 ' 5 M decreases 
rapidly, and a fraction of these galaxies enters into a qui- 
escent phase (massive galaxies with low SSFR in Figure 
6 and 7). As a result, the contribution of these galaxies 
to the SFRD becomes relatively small by z ~ 2. At 
z < 2, the decrease of the SFR in each galaxy with 
Mstar = 10 10-11 M overwhelms the increase in the num- 
ber density, which drives the overall decrease of the cos- 
mic SFRD with time. 

In Figure 12 and 13, we show the contribution of 
galaxies with different SSFRs to the cosmic SFRD in 



Figure 13. The same as Figure 12 but for the case where 
SFR correct ed-uv is used for all the sample. 

the cases with SFRjr + uv an d SFR correcte d-uv? respec- 
tively. It is seen that in galaxies with M s tar = 10 10-11,5 
M , those with a relatively narrow range of SSFR (< 1 
dex) dominate the cosmic SFRD, especially in the case 
with SFRir + tjv- Furthermore, the SSFR of galaxies 
which dominate the SFRD systematically increases with 
redshift. In this context, major star formation in the 
universe at higher redshift seems to be associated with 
a more rapid growth of stellar mass of galaxies. At 
2.5 < z < 3.5, in particular, the high-SSFR popula- 
tion with SFR/Mstar ~ 10 Gyr -1 , which are relatively 
small in number (Figure 8), dominate the SFRD. If these 
objects would be the 'starburst' galaxies with distinc- 
tively higher star formation efficiency than ordinary disk 
galaxies as discussed in the previous section, these galax- 
ies are expected to grow their stellar mass rapidly and 
could cause the increase of the number density of mas- 
sive galaxies at 1 < z < 3. There might be a transition 
in the dominant component of the star formation in the 
universe between z ~ 3 and z ~ 1; from starburst in rel- 
atively small number of galaxies at z ~ 3, which drives a 
rapid growth of stellar mass, to star formation in many 
disk galaxies at z ~ 1, which leads to a more gradual 
growth of stellar mass. 

7. SUMMARY 

We studied SFR as a function of M s tar for galaxies 
at 0.5 < z < 3.5 in the GOODS-North field, using the 
if -selected sample from the MODS. The very deep NIR 
data of the MODS allow us to construct the stellar mass 
limited sample down to io 9 - 5-10 M even at z ~ 3. The 
rest-frame 2800 A luminosity and the MIPS 24 fim flux 
were used to estimate the SFRs of the sample galaxies. 
Using these SFR indicators, we showed the two cases of 
the results, namely, one case where SFRir + tjv is used 
if a galaxy is detected in the MIPS 24/im image and 
SFR CO rrected— uv is used for those without the 24/im de- 
tection, and the other case where SFR correcte d-uv is used 
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for all the sample. 
Main results in this study are as follows. 

• More massive galaxies tend to have higher SFRs 
at M star < 10 10 - 5 M in all redshift ranges. At 
M s tar ^ 10 10 - 5 M , there are galaxies with rel- 
atively low SFRs similar with low-mass galaxies, 
while galaxies with higher SFRs which follow the 
trend at lower M sta r also exist. 

• The median and average SSFRs of galaxies is nearly 
independent of M sta r at M sta r < 10 10 5 M at z > 
1, while they slightly decrease with M sta r at 0.5 < 
z < 1.0. At M star > 10 10 - 5 M , the median and 
average SSFRs decrease with M star . 

• The SFR (SSFR) distribution at a fixed M star shifts 
to higher values with increasing redshift. More 
massive galaxies show the stronger evolution of 
SFR at z > 1. 

• At 2.5 < z < 3.5, the bimodality in the SFR and 
SSFR distributions at each mass is seen. The bi- 
modality is the more simple form in SSFR and can 
be divided into the two populations by a constant 
SSFR of ~ 2 Gyr -1 . One population consists of 
galaxies with SFR/M star ~ 0.5-1 Gyr -1 , and the 
other population is the high-SSFR population, to 
which galaxies with SFR/M star ~ 10 Gyr -1 belong. 
These high- and low-SSFR populations might have 
different modes of star formation, namely, star- 
burst and more continuous star formation in disks. 

• Galaxies with M star = 10 10-11 M dominate the 
cosmic SFRD at 0.5 < z < 3.5. The evolution 
of these galaxies seems to drive the evolution of 
the cosmic SFRD, especially at z < 2.5. The con- 
tribution of galaxies with M star = 10 11-11 - 5 M 
increases strongly with redshift at z > 1.5 and be- 
comes significant at z ~ 3. The contribution of 
galaxies with M star < 10 10 M is relatively small 
at any redshift. 

• In galaxies with M star = 10 10-11 - 5 M star , those with 
a relatively narrow range of SSFR (< 1 dex) dom- 
inate the cosmic SFR at 0.5 < z < 3.5 especially 
in the case with SFRtr+tjv- The SSFR of galax- 
ies which dominate the SFRD systematically in- 
creases with redshift. At 2.5 < z < 3.5, the high- 
SSFR population dominates the cosmic SFRD. Ma- 
jor star formation in the universe at higher redshift 
seems to be associated with a more rapid growth 
of stellar mass of galaxies. 

Systematic errors in the SFRs estimated from the 24 
fim flux and the optical-NIR broad-band SED (Section 
4.3) could affect some of the above results. More pre- 
cise measurements of SFR are essential to conclusively 
confirm the results in this study and to reveal the stellar 
mass assembly histories of galaxies at 1 < z < 3. For ex- 
ample, NIR spectroscopic surveys with sufficiently large 
sample for the stacking analysis and/or long exposure 
time enable to estimate directly the dust extinction for 
the nebular emission lines from the Balmer decrement, 
and to determine the dust-corrected Ha luminosity of 
star-forming galaxies at that epoch with high accuracy. 



The forthcoming Herschel and ALMA observations will 
also provide strong constraints on the shape of the far-IR 
SED of these galaxies, which allows us to determine the 
total IR luminosity more precisely. 
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